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ABSTRACT 


A comparative study is made on the realizability conditions of 
various special cases as well as the general case of the positive-real 
biquadratic functions. Computer programs are written for this study to 
illustrate the realizability regions of zeros (or poles) of the functions. 
A limited study is also made on the sensitivity of the function to 


variation of element value in a realization. 
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I. INTRODUCTION 


The problem of realizability of biquadratic functions as driving 
point immitance has been studied by a number of investigators [1] - [11], 
with realizability conditions given for various special cases as well as 
for the general case of biquadratic functions. It has been proved [10] 
that, given a pair of poles, the location of the zeros are limited by 
two curves, one of which is the locus of the zeros such that the biquad- 
ratic function is a minimum positive real function. 

In Chapter II of this paper, various realizability conditions are 
reviewed and illustrated. FORTRAN programs are written in Chapter III 
for illustration as well as comparison of different cases. 

Many forms of realizations of positive-real biquadratic functions 
are well known. Some of these forms are utilized, in Chapter IV, to 
study the sensitivity problem related to the realizability of the 
function. The method of graphical representation by Chan [13] is applied 
in this study. A computer program is written for the illustration of 
the variation of sensitivity phasors [13] as functions of element values. 

Finally, in Chapter V, a discussion on the results is made together 


with some suggestions for further studies. 





II. REALIZABILITY CONDITIONS 


A. INTRODUCTION 


Consider a general biquadratic function 


Z 2 
St 20's + 
Z(s) = a = eisai (Stee (1) 
re (s + p')(s + p") 
P P 
or 
ae + long + ne (oe - me) + jw2o 
: _ z Zz 2 2 
2(ju) = —>——2——2 - —2_, ; (2) 
-w + 20 jw + w (w' - w ) + jw2o 
P Pp P P 


It has been investigated by Chan and Phung [10], that a biquadratic 
function Z(s) as given in equation (1) with conjugate-complex zeros in 
the open half plane is realizable as a driving point immitance with 
passive one-port if and only if the zeros are located either on the 
curve 2 or Cc, or inside the area delimited by curves ee and C5, whose 
polar equations are given respectively by 


t ? s 
max w= Ee + en cos uy + ae cos Dake cos a + “» | (3-a) 


and 
? 


' ? 4 
Minw = ue + a cos oe - 2 E cos oe cos Oe =e > | (3-b) 


' ? 
where oe denotes the compliment of the argument of the zero z , such 


that ’ 


So, =. cos > . 
Zz z Zz 








Figure l 


B. NECESSARY AND SUFFICIENT CONDITIONS 
To obtain the real part of Z(jw) of equation (2), multiply and 


divide this equation with the conjugate of the denominator, so that 


Ga - 0") (we - ne) + 4o ou 
Re Eo = cal ana iel ie aaa a : (5) 
on -w ) + te 


Z(s) is positive real if the numerator of expression (5) is nonnegative 


for all values of w, that is 


nie) a i - Poe ~ a - boo) + ote 
Z Pp zp 


oes (6) 


Now we have the following two conditions [10]: 
Condition-1 

N ee has complex or equal positive real roots. 
Condition-2 

N (a has nonpositive real roots. 


; ee aa ; 
Condition-1 implies the discriminant of N(w') is negative or zero 


D= ie ~ ie - 40 a - - nie 20 (7) 
Zz Dp 20 Zp 





or after factoring we get 


2 2 
le, ee) - io, | lw, ~ ay - 2%| <0; (7-a) 


thus condition-l is satisfied if 


(w, - ws)” $ too, < (wu, tw) , (7-b) 

The right hand inequality is shown as strict since inequality corresponds 
to real roots which would in this case be nonpositive, which is a 
violation of condition-l. 

Next, condition-2 implies the following two conditions simultaneously: 
Condition-2a 

Dy,.0O. 

Condition-2b 


Oe ~ w* —t4coo< 0 . 
zZ Pp 2p 


It can be seen that both conditions 2a and 2b are satisfied if 


Pe 
seo, 3 (w, + a) : (8) 


With equations (7-b) and (8) we obtain the necessary condition for positive 


realness of Z(s) to be 


Qe < icc. (9) 
Zz Pp Zap 


Now suppose that equation (9) is satisfied with the equality sign. 
; 2 
Z(s) is then a minimum positive-real function since N(w ) has a double 


positive root Wis where 


Ww. = WW : (10) 
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Next suppose equation (9) is satisfied but not necessary as an equality, 


then 


(ow - 4) 
oO, 2 mae : (11) 


The necessary and sufficient condition obtained in the above discuss- 
ion can be conveniently and geometrically interpreted. 
Suppose that the zeros z' and z" are conjungate complex with negative 


real parts and the poles p' and p" are similarily located in the open 


left half plane. 


Let 
' 


Oo, = w, cos >) (12-a) 


' 
where pe denotes the compliment of the argument of the zero z’. 


If w 2w_, then we let w =w tu (12-b) 
iz P Zi 


(12-c) 


i) 
Ee 
1 
G 


If w <w _, then we let wo 
Zz Pp Z 


substituting equations (12-a) and (12-b) into equation (9), and after 
rearranging we get 


2 ! ' 
us 7 oe (cos oe) i ace cos >, < 0 (13-a) 


and substituting equations (12-a) and (12-c) into equation (9) we get 


' 


2 ' 
- : 13=b 
u, + “en (cos >.) u ona cos o>. < 0 (13-b) 


Both polynomials in the left hand sides of the inequalities (13-a) and 


(13-b) have a positive real root and a negative real root, so that these 


lat: 





two conditions will be satisfied if uy and u, are nonnegative and smaller 
than or equal to the positive-real root of the left hand side polynomials 
in (13-a) and (13-b) respectively. 

Thus solving for the positive real root and then substituting the 
result into equations (12-b) and (12-c) gives, for each value of _ 06.4 
the maximum and minimum values of Ws for which the necessary and suffi- 
cient condition for the positive realness of Z(s) is satisfied. 


Alternatively, we substitute equations (12-b) and (12-c) into equation 


(9) and get 
2 
u- mee, <0 
and 
ie =iog Go< O . 
i Zap 


Solving these two equations for the positive real root and substituting 
the results into the expression for Wo» gives the maximum and minimum 
values of w for which the positive realness of z(s) is satisfied. 


Thus we have in this case the two equations 


maxw =w +2 vo o (14-a) 
z P Zz DP 
and 
minw =w -2 y¥0o0 ; (14-b) 
z Pp Zp 


Note that these equations will give the same results as equations (13-a) 
and (13-b) when poles of Z(s) are complex, and will also give the 


realizable area of zeros of Z(s) when the poles are real. 


Thus given Z(s) with complex conjugate poles in the left half plane, 
or with negative real poles, Z(s) will be realizable if the zeros of Z(s) 


are complex and lie on or within the curves of equations (14-a) and (14-b). 
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In the special case where Z(s) has real poles and complex zeros, 
it has been shown that z(s) can be realized by a five element series 
parallel circuit [11]. 


The key restriction in this case is given by 
p" (Gu ps p') =, Cz _ p”) Cz. = a) > 0 C15) 


where p', p", z' and z" are poles and zeros of 


(s + z')(s + 2") 


2(8) = (s+ p(s + PY) 


ana p' > p'. 
If we substitute o's and w's for the p's and z's of (15), then we 


get 


z 

we also obtain 

20 o0 + 20 Ye = a = te - ve > O 
zp z p Pp 


rearranging terms, we have 


oo a a < 20 6 +2 oe = i. 
Zz Pp zp Z Pp p 


Now since 


13 





then we get 


(Ww -w i. Sc Oe 620) Vea - na - 20 w _ 
- Pp Zp Z p P zp 


Comparing this condition to that for the general case in equation (9), 


we obtain 


Zona 20 Jo” = na < 40.00 ; 
Zp Z Pp P zp 


and since ae a 2 and are all positive numbers, the result becomes 


(ow, - ws)” < 2050. + 20, Yo = WF = ae < ae 
so that this is indeed more restrictive than the general case and we 
should not be surprised to find the area of realizability for this case 
enclosed by that of the general case. 

In fact this argument could be used as a proof of the realizability 
of Z(s) given in (15), since the area of this case is always enclosed 
by that of the general case as shown above. 

The construction of the realizability region in this case will be in 


accordance with reference [11], that is 
2 3 
po M=re= p> (p" a DD ; 
where p" > p'. 


C. SAMPLE ILLUSTRATIONS 


Compl pol =O 
(a) omplex poles Co, > 


Let oF = 3 
Pp 


and w = 5 
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Then from equations (14-a) and (14-b) we get 


Ww max =w + 2Vv¥0 a 
P zp 


5+ 2v30_ =5+ v120, 3 


and 


W Min =w - 2vV0 ao 
z Pp zp 


z 
We obtain the limitation of ce as follows: 
oS < uw for complex zeros. 


If we let wu max = Sits then we get 


S 
0 


5+ vi2o5 


20.8 = — maximum value for os for w, max. 


Similarly if aa nin = Oo: then 


oon 2 = O8 maximum value for J for a min. 


(b) Real Poles 


G7 2.55 and 
wy = 2 with poles at (-l, -4). 


Then we obtain 


WO max 
Zz 


we +r 2¥ a0 
Pp Zp 


2 + Y100 
Zz 


se [lal aed 





and 


€ 
tt 


2 + Y100 
Zz 


Lse7 = o> max o, EOF w max and complex zeros. 


From the examples given in this section it can be seen that the 
realizability region of a biquadratic function is dependent upon the 
restriction we place on the function. That is, as we place more 
restrictions on the function and proceed from a general function with 
complex poles to one with real poles which is realizable by certain 
technique, we see that at each step the realizability region decreases. 
This decrease of freedom from the additional restriction is to be 


expected. 
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IIT. COMPUTER STUDY ON REALIZABILITY 


For a given location of the complex pole pair i and a, specified) 
in the open left half plane, C, and Cc, are simply plots of w versus ,. 
Thus the computer program was written to accept o, and w_ as inputs and 
to generate C. and C; on either the line printer or the Calcomp plotter. 
Some sample outputs are shown in Figure A-1.A and A-1.B in Appendix A. 


The curves are for a normilized set of pole locations with 


and 
ba = é < 90°, 
P 
' ' 
where . is the compliment of the argument of the pole p . 

Varying 5 from 0° to 90° corresponds to shifting the poles along 
a circular path from the negative real axis to the jw axis. 

The program was then modified to produce a set of curves Cy and Cc, 
for various pole locations. Given a specific oie of poles, two types 
of plots may be obtained, depending upon the requirements of the problem. 

If an accurate plot of the realizability region is required, the 
program can be used so that a plot of the realizability region is 
required, the program can be used so that a plot similar to that shown 
in Figure A~1.A is obtained. Otherwise if a quick estimate is all 
that is desired, a plot similar to the one in Figure A-1.B can be used. 

Other outputs were obtained for various pole locations. 

Results obtained for the following cases are included, with 


program listing, in Appendix A and B respectively: 


1. Complex poles with oe moves toward a constant oe 


17 
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7 


Se 


Complex poles with 2 moves toward a constant Soe 

Complex poles that are moving toward the origin. 

Real negative poles with one of which moving toward the origin 
while the other remains stationary. 

Real negative poles with one of which moving toward infinity 
while the other remains stationary. 

Double real poles moving from the neighborhood of the origin 
toward infinity. 

Pure imaginary poles moving toward the origin. 


Pure imaginary poles moving away from the origin. 


The above cases are now described in the following sections. In 


each case a figure is included showing the pole locations with their 


values shown in the accompanying table, in which the location of the 


corresponding computer output plots are indicated. 
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Case l 
This case illustrates the zeros region for given complex poles with 


o moves toward a constant w. 
P P 





Figure 2. 


Five different poles ( conjugate poles are not shown ) are given in 
accordance to show the zero regions for each Oo, as My remains constant. 


The locations of each pole are tabulated in Table l. 


Table 1 
Pole oO w Computer Plot 
ee 2 pee eee a il pape pe SS 
Cc, 0.999 1.0 Figure A-1.A,B 
C, 02707 ie Figure A-2 
Ce Oc1 1.0 Figure A-5 
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G@ase 2 





Biguniens >. 


This case illustrates the zeros region for given complex poles with 
Wp Moves toward a constant og. 
P 


Five different pole locations are shown in Table 2. 


Table 2 
Pole 0 Ww Computer Plot 
P P 
a 1.0 1.05 All plots are 
CS eG 2.0 shown in 
C 1.0 230 
4 
C 1.0 8.0 
5 
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Case 3 





Figure 4. 


This case illustrates the zeros region for given complex poles that 
are moving toward the origin. Five different pole locations are 


shown in Table 3. 


hapless 
Pole ¢., ») Computer Plot 
P 
Cy 4.0 30.652 All plots are 
C, 2.0 15.326 shown in 
C, Onl 0.766 
Cy q. 01 0.077 


Zi 





Case 4 
This case illustrates the zeros region for given real negative poles 
with one of which moving toward the origin while the other remains sta- 


tionary. 





Five different pole locations are shown in Table 4. 


Table 4 

ee cE iE 

Condition Pole PD ole 2 C, Ww, 
ci 3.0, 0 4.0, 0 6.0 52057 
C, 60), 0 InOy, aS 2.020 
C, 6.0, 0 G25;,..0 4.25 2.50 
C, 5.0. 0 Onl, 4.05 1.414 
C, o.0, U0 0..05;0 4.025 0.632 


All plots are shown in computer output Figure A-9. 





Case 5 
This case illustrates the zeros region for given real negative poles 
with one of which moving toward infinity while the other remains stati- 


onary. 





Five different pole locations are shown in Table 5. 


Table 5 


Condition Pole PL Pole ES o wy 
e 8.0,0 2.0, 0 5.0 4.0 
S 15.0,0 2.0, 0 one) 5.477 
S 30.0,0 20s 16.0 "7.746 
Tp 50.0,0 2.0, 0 26.0 10.0 
c. 100.0,0 Ze Ona 0) 56.0 14.142 


All plots are shown in computer output Figure A-10. 


Z3 





Case 6 


(jw) 





Figure 7. 


This case illustrates the zeros region for given double real poles 
moving from the neighbourhood of the origin toward infinity. Five 


different pole locations are shown in Table 6. 


Table 6 


Condition Pole P,=P. ae a, Computer Plot 
Cc) 1.0, 0 a0 1.0 AIl plots are 
Co 200. 0 2.0 22.0 shown in compu- 
C, 5050 0 2120 5.0 ter output 
Cy, 10.0, O 10.0 10.0 
C. 25.0. 0 a) 25.0 
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Case 7 


(jw) 


Figure &. 


This case illustrates the zeros region for given pure imaginary poles 


moving toward the origin. Five different pole locations are shown in 
Table 7. 
Table 7 
Condition ae 2. Computer Plot 

of 0.0 3.0 All plots are 

C, 0.0 1.0 shown in 

C, 0.0 0.5 Figure A-12. 

C, 0.0 0.1 

C, 0.0 0.05 


Z5 





Case 8 





Figure 9. 


This case illustrates the zeros region for given pure imaginary poles 


moving away from the origin. 


Table 8. 


Table 8 


Five different pole locations are shown in 


Computer Plot 


Condition CG 
P 
C) 0.0 
Cy 0.0 
C, 9) 
C, 0.0 
Ce 0.0 


10.0 


fo 
q) 
° 

© 


100 .0 


200.0 


26 


All plots are 
shown in com- 
puter output 


Figure A-13. 








A. INTRODUCTION 


IV. SENSITIVITY STUDY 


It has been shown [12], that in a linear and not necessarily 


reciprocal network, with only one independent source and initially at 


rest, the ratio of the Laplace transform of the current in any branch 


or of the voltage across any node pair to the Laplace transform of 


either an independent voltage or current source is given by 


Tx) 


where 


T(x) 


T (0) 


T (=) 


ut 


il 


W T (0) ters Eke) 


Wa x (16) 

an adjustable immitance. 
the network function relating a response to an 
exitation. 
Lin T(x) = the network function evaluated with 

x- 0 
x= 0, T (0) is assumed to be finite. 
Lim T(x) = the network function evaluated when 


>_< ———- © 


x = 2%, T (2) is assumed to be finite. 
the Thevenin immitance (independent of x) seen looking 
back into network from the terminals of the adjustable 


immitance. W is assumed to be finite. 


Equation (16) is generic in form in that x can be either an impedance 


or admittance. 


point immitance. 


In this paper the symbol T(x) is used to denote a driving 


ay 





B. GRAPHICAL REPRESENTATION OF T(x) 
Consider equation (16) for the case of a driving point immitance 
function as shown in Figure 10. 


Dividing numerator and denominator of equation (16) by W, we get 


x 
70) + —~ T (=) 


T(x) = iM 


x 
eee 
W 


which can be rewritten as, 


x 
Ree) TO) TO) 
T(x) = 
x 
1+ — 
W 
or 
tet (2) 
Tene) (17) 
x 
1 +— 
W 
Let 
Ueyta(e) + V (18) 
where 
te(O)) = T=) 
ee ss ; (19) 
x 
1 +— 
W 


28 








BRagure 10. 


Driving point immitance with variable 
component x. 





Figure ll. 


Graphical representation of T(x) phasor. 
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From equation (19) we have 


x 
VEV mee = T (0) = =) (20) 


we note that T(0) and T () are fixed quantities. 


Now treat each quantity in (20) as phasor, i.e. 


<j 
i 


lv| exp. 58, 


and 


= 
u 


\w| exp. j8, 


then the left hand side of equation (20) represents the sum of two phasors. 
If x is real, then the phasor 
x V 
wa OX. |——| exp. j(6, - 6.) 
W W 1 Z 


is lagging behind the phasor V by an angle 


8 = 85 = 84: 


Since T (0) and Ee) are two fixed quantities, the right hand side 
of equation (20) is constant so that the left hand side of (20) is 
constant. This requires that the angle must be constant (as x varies). 

This implies that the tip of phasor V describes the arc of the 
circle, of which the quantity T (0) ~ T (eo) is a chord as shown in 
Figure ll. 

To construct the circle we proceed as follows [13]: 

ie Plot Eee) and T (0) as two phasors from the origin. 

2. Drawa straight line D, through the tips of the two phasors in 


step l. 


30 





3. Calculate the angle and draw a straight line dD, from the tip of 
T,.(0) forming angle with Dy having the value -@ measuring clock- 
wise from dD): 
Therefore dD, is tangent to the circular locus of T,- 
The centre of this circle, 0', is located at the intersection of 
the two perpendicular lines erected on the chord in the middle of 
T (0) - ae) and Dy at the tip of T (0), as shown in Figure 12. 


Knowing the locus of T(x) as x varies from zero to infinity, we 


can graphically determine T(x) at any value of x. 


C. GRAPHICAL REPRESENTATION OF SENSITIVITY PHASOR 


If we define the sensitivity as 


Hoe. d T(x) 
S = 
x dx 


then with equation (16) we can obtain, [13] 


T (0) - T (=) 


my 
= - 2 
x 
will+-s- 
W 
or 
x gi (x) 
ee 
x T (0) - T (2) 
vV-—- 
W 
where 


T (0) - T.(@) 
x 
oo 
W 


Vv = 


31 





cus of T(x) 









RIguGbew 2. 


Graphical representation of locus of T(x). 





Brgure 13. 


Graphical representation of sensitivity 


phasor S. 
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It can be shown [13] that the sensitivity of a function may be 


graphically represented by a phasor such as PS as shown in Pigure 15. 


D. SAMPLE ILLUSTRATION 

As an illustration of the above discussion in this chapter, let us 
consider a biquadratic p.r. function as an RC driving point admittance, 
which can be realized in Foster 2 form as shown in Figure 14. 


Let 


fife t 3 — Gener 


Y(s) = 5 
s + 6s + 8 


which can be realized in Foster 2 form with the element value as given 
in Figure 5, in which all resistance and capacitances are measured in 
Ohms and Farads respectively. 


Now .let Ry change by the amount AR, such that 


0 < AR < @ 
Thus, Figure 15 is modified as shown in Fgiure 16. 


We can calculate T, (0), Ee) and W as follows: 


1 1 1 
Po 4 9@ Se 
x Ry Ry + 1/sC, Ry + 1/sC, 
125 09375 s Soe ee 
eter 5aen 2, ,, © 2 
5 ‘ s + 6s + 8 
aS reel Zz 
eee = Soe, * a i/sc 
oO i a) 
: oma es a 3 (25) 
ewer " 45 +3 
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Figure 14. 


The realization of Foster-2 form of 
driving point admittance Y(s). 


Rl =4 R 


te 
“- 
nw 

eel 


Pigure 25. 


The realization of biquadratic p.r. driving point 


admittance in Foster 2 form, where 
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oy pet a * ¥(s, AR==) 
_ 10.667 4s+ 8 
= 2.667 + 754 “ ae 
14.5 ao + 40 s + 99 
—- 7 (26) 


758s +95 


From expression (26) we obtain the expression —— to be 


2 
* ms mie 5 a x+9sx (27) 
14.5 s + 40 s + 99 


Substitute equations (24), (25), (26) and (27) into equation (22) we get, 


Bee ee Siete 


2 
gt (*) eS 0 i ok | 2: » (28) 


NARS e> oeOne 2-90 E eh xcs 0 xenale 
Tee one 14.5 s¢ + 40 s + 99 


then for simplification we write equation (28) above as, 


T(x) 
S = eee 
x 


N(s,x) 


eS) (2?) 


where N9s,x) and D(s,x) represent the numerator and denominator of the 
equation (28) respectively. 
In the frequency domain, the expression in equation (29) can be 


written as, 


() j 
eGo) = ee (30) 
x 


a5 





row 


! ; 4R 
| (=> 
Ry : 

Y(s)! 

| R 
| Oo R 
\ 2 
Cy C2 
| 
Figure 16. 


Driving point admittance Y(s) with variable resis- 


tance 4R. 


'W 
J {| ----- = 
(i ale ARI 
| | 
| | } 
! a eal ee 
7 | R | 
oO 
t | | Ro: 
1 ; geen ote ae ) 
| Cy C. 
| 
Figure l7. 
Driving point admittance Y(s) with variable resis- 
ie 
tance Rs = R, try ot 6 
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Re {N(jw,x)] + j Im [N(jw,x) ] 


(30a) 
Re [D(jw,x)] + j Im [D(jw,x) ] 


where Re [N(jw,x)], Re [D(jw,x)], Im [N(jw,x)] and Im [D(jw,x)] are 
real and imaginary parts of equation (30) repectively. 
From equation (30a) we can obtain the magnitude of that sensitivity 


function. If we write equation (30a) in exponential form we have 


TGs) = Sy LG) 
Sy Gusx) = [S| exp (ty - 45) (31) 
where 
Leg) _ : : 
ae = magnitude of equation (30a) 
{ : 14 : Z 3 
Re [N(jw,x)]}" + {Im[N(jw,x)]} 
: 2 : Z 
{Re [D(jw,x)] ~ + {Im[D(jo,x) ]} 
oo = the argument of the numerator of equation (30a). 
$y = the argument of the demominator of equation (30a). 


Be COMPUTER STUDY ON SENSITIVITY 
In this section, the sensitivity of the function T(x) with respect to 
the variation of one of its components will be investigated by computer. 
Consider the general expression Y(x) in equation (16), which is repeated 
here for convenience. 
W T,(s,x =0) +x Ts; oe Sa) 


T(x) = —“———_____*—_____.. (16) 
Wt+s 


The sensitivity function of this equation is 


x 7 x 2 
W [2 += ] 
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S (22) 





As an example let us consider the function Y(s), as we derived in section C, 


(Sees 
(s + 2)(s + 4) 


Y(s) 
and the realization of this function in Foster-2 form is shown in Figure 
ay 
For our investigation, let Ry change with the amount of AR in series 


with R, and let 


' = 
R, R, po AR. 


q 
To obtain the sensitivity sie a ) 
Z 


we use the same calculation from equation (30). 


in terms of frequency and Ro! 


The computer output in Appendix B shows us the plot of the magnitude 
of sensitivity function (30) versus R,'- For this purpose we use equation 
(31) as reference. 

Various cases are investigated using computer, resulting in plots 


of the magnitude of equation (31) versus R,' for different frequencies. 


Z 


38 





V. CONCLUSION 


In this paper, a comparative study has been made on the realizability 
conditions of a general biquadratic function. The graphical method of 
Chan and Phung [10] was used for this study and computer programs are 
written for varifications. General observations can be made as follows: 

1. The realizability for the zeros expands considerably as the poles 
move from the immaginary axis to the real axis. This results both from 
C, moving away from the origin and Cc; moving toward the origin. 

2. The limits on C, and C, are found easily from equations (3a) and 
(3b). As a approaches zero, corresponding to the poles moving toward 
the jw axis, both C, and C; tend toward the mE cie w= . for alu oo: 
As can be seen from Figure A-6, C, and Cc, separate very rapidly for 
small movements of the poles away from the jw axis. We can conclude that 
when oe = 0 we require w, = oe in order for Z(s) to be realizable with 
complex zeros. This can also be seen easily from equation (9). 

3. As c, approaches Boe corresponding to the poles moving toward 
the real axis, Co approaches a value of a. on the real axis, and 
Cc, similarly approaches a value of eee as we can see in Figure A-1.A. 

A brief study was made on sensitivity functions [13], some computer 
results are obtained. However, due to time limitation, a complete 
investigation in sensitivity was not possible and is left as a 


suggested topic for further studies. 
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APPENDIX A 


COMPUTER OUTPUT 


A. ZERO'S REGION 

The illustration of zeros regions will be shown in this appendix. 

The computer outputs show us the regions for various given omega 
and sigma poles. 

Various cases are explained as follows: 

Figure A-la and A-lb show the zeros regions for given sigma pole 
= 0.999 and omega pole = 1.0 using Calcomp and PLOTP in the computer 
program respectively. 

All the following plots show us the computer outputs by using Calcom 
in the computer progran. 

Figure A-2, A-3, A-4 and A-5 show us the plots of zeros regions for 
given sigma poles equal to 0.707, 0.55, 0.25 and 0.1 respectively with 
the amount of omega pole equal to 1.0 each. 

Figure A-6 shows all plots given in Figure A-la, A-2, A-3, A-4 and 
A-5 in one computer output. 

Figure A-7 shows the plots of zeros regions as pole moves at sigma 
pole constant. 

Figure A-8 shows the zeros regions as pole moves to the origin. 

Figure A-9 and A-10 show the plots of zeros regions as pole moves 
to the origin and to infinity respectively. 

Figure A-1l shows the plot of zeros regions of a double real negative 
poles as both move along the negative axis. 

Figure A-12 and A-13 show us the zeros region plots of a given pair 
of pure negative poles, as both move to the origin and to the infinity 


respectively. 
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B. SENSITIVITY FUNCTION 

Various cases of sensitivity as a function of Ry are illustrated in 
the computer outputs for different particular frequencies. 

The computer outputs show us the plots of sensitivity versus R'4 
and the tabulations of this plot. 

The resistance R', is select to be varied from the amount of 1 Ohm 
to 50 Ohms with increment to 1 Ohm. And the frequencies are choosen 
beginning at 1 rad/sec. up to 10,000 rad/sec. 

Figure B-l1, B-2, ...... B-9 show us the sensitivity functions as 
R', moves from 1 to 50 Ohms, at the frequencies of 1, 5, 10, 50, 100, 
500, 1000, 5000 and 10000 rad/sec. The tabulations of those plots 
are also shown in Table B-1, B-2 ....... B-9 respectively. 


Figure B-10 shows us five different plots of sensitivity functions 


illustrated in one computer output. 
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APPENDIX B 


COMPUTER PROGRAM 


A. ZEROS REGION 
This program will plot the realizability region of general bi- 
quadratic function of the form 


= + 2*SIGMAZ%s + OMEGAZ? 





Z(s) = 


eo + 2*SIGMAP%s + OMEGAP* 


The input to this program is via punched card. One card is required 
for each plot. The input variables are Sigmap and Omegap specified in 
columns 1 - 20 as two fields of F 10.3 format. 

The output is a plot of allowable region for zeros of Z(s) for a given 
pair of poles. The criterion for the plot is that Z(s) must be positive 
real. 

The output may be either one (or both) of two types of plots depend- 
ing on the needs of the programmers. If a fast solution is desired the 
"CALL DRAW" statements must be removed to enable the program to be run 
on QUICKRUN. The result will be a rough plot on the line printer 
output. 

If maximum accuracy is desired a CALCOMP plot may be obtained by 
making a regular FORTCLG run with the "CALL DRAW" cards in the deck. 

In this case a line printer plot will also be generated, unless the 
"CALL PLOTP" cards are moved. For the CALCOMP run, "REGION.GO-100K" 
must be specified on the EXEC card. 

The routines DRAW and PLOTP are called by this program. Both are 
in the MPSLIB at NPS. Neither a copy of the decks nor special job 
control language is needed to run this program on the NPS IBM 360. 
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of general biquadratic function. 





The computer program is shown in Figure A-14. 


B. SENSITIVITY FUNCTION 
This program will plot the sensitivity of general biquadratic 
function with respect to the changing of its components at a particular 
frequency. 
The equation used in this program is, 


W (T(jw,R, = 0) - T(ju,R} = INF.) 


T(jw,R 2 
SR 


2) (ju,R,) ee. 2 
(W + Ro) 


where T(jw,R,), W, T({w,R, = 0), T(jw,R, = INF.) and Ry correspond to 
Figure 17 in Chapter IV. E. 


The output is a plot of the magnitude of the sensitivity function 
above versus Ry. 

The input of this program is via punched card. One card is required 
for each plot. The input variable is frequency in rad/sec. specified 
in columns 1 - 10 as one field of F 10.3 format. 

The output is a rough plot on line printer output. 


The routine "CALL PLOTP" is called by this program. 


The computer program is shown in Figure B-1ll. 
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FIG B-11 


Computer program for sensitivity function of 


given particular biquadratic function. 
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